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Abstract 

We consider a stochastic delay differential equation driven by a general 
Levy process. Both, the drift and the noise term may depend on the past, 
but only the drift term is assumed to be linear. We show that the segment 
process is eventually Feller, but in general not eventually strong Feller on 
the Skorokhod space. The existence of an invariant measure is shown by 
proving tightness of the segments using semimartingale characteristics and 
the Krylov-Bogoliubov method. A counterexample shows that the stationary 
solution in completely general situations may not be unique, but in more 
specific cases uniqueness is established. 

1 Introduction 

Stochastic delay differential equations, also known as stochastic functional dif- 
ferential equations, are a natural generalisation of stochastic ordinary differential 
equations by allowing the coefficients to depend on values in the past. When 
only the drift coefficient depends on the past, main stochastic tools and results 
for stochastic ordinary differential equations can be applied, for example by re- 
moving the drift via a change of measure. If the stochastic perturbation depends 
on the past, however, surprising new phenomena emerge, see Mohammed and 
Scheutzow [22] for a discussion on flow and stability properties. 

*0. van Gaans acknowledges the financial support provided through the European Commu- 
nity's Human Potential Programme under contracts HPRN-CT-2000-00100, DYNSTOCH and 
HPRN-CT-2002-00281. 
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The main purpose of the present work is to investigate the stationarity of delay 
differential equations driven by Levy processes of the form 



L denotes a general Levy process, the drift term is obtained by integrating past 
values with respect to a signed measure \i and the nonlinear coefficient F(X) 
depends on (X(s) : s £ [t — a,t]) at time t, see Section 2 for details. We do 
not consider a nonlinearity in the drift in order to concentrate on the effects 
of the nonlinear noise term, which is facilitated by a variation of constants for- 
mula. While the solution processes are not Markovian anymore, one can retrieve 
the Markov property by regarding segments of the trajectories as processes in 
a function space. The delayed noise term causes a fundamental degeneration of 
the segment process: we show that the Markov semigroup is not Feller and not 
eventually strong Feller, but eventually Feller. Consequently the uniqueness of 
an invariant measure can not be derived by the strong Feller property. 
Stationarity results for Levy-driven stochastic differential equations, even in 
the non-delay case, are not so widespread. Non-Gaussian stationary Ornstein- 
Uhlenbeck processes have been attracting increasing attention recently due to 
their use in financial modelling and the relationship with self-decomposable dis- 
tributions, cf. Barndorff-Nielsen and Shephard [4]. Invariant measures for ordi- 
nary differential equations with a nonlinear drift term and additive stable noise 
have been studied analytically by Albeverio, Rudiger and Wu [1], but our general 
results, even when specified to the non-delayed case, seem to be new. 
The question of the existence of stationary solutions of stochastic equations with 
delay goes back to the 60s in the work of Ito and Nisio [14]. They have proved the 
existence, but not the uniqueness of a stationary solution for Wiener-driven delay 
differential equations under the condition that the drift is obtained by a delayed 
perturbation of a stable instantaneous feedback. For a more general non-linear 
drift functional and additive white noise, Scheutzow [28] derived sufficient condi- 
tions for the existence of an invariant probability measure in terms of Lyapunov 
functionals. For a similar approach and connections to stochastic partial differ- 
ential equations see Bakhtin and Mattingly [3]. For the Levy-driven equation 
(1.1) with constant F, that is additive noise, Gushchin and Kiichler [13] have 
established the existence and uniqueness of stationary solutions. 
Our work is based on analyzing the segment process in a function space and it is 
therefore closely related to results for stochastic evolution equations in infinite- 
dimensional spaces. In case of additive noise an extensive literature for the sta- 
tionarity of solutions of stochastic evolution equations exists, see Da Prato and 
Zabczyk [10]. Much less is known for non-additive noise, see for example Chow 
and Khasminskii [8] for some general results. An infinite-dimensional analogue 
of equation (1.1) driven by a Wiener process is considered by Bonaccorsi and 
Tessitore [6]. They obtain a stationarity result for small Lipschitz constants by 
a fixed point argument. 

To prove the existence of a stationary solution of (1.1) under rather general 
conditions, we consider the segment process with values in the Skorokhod space 




(1.1) 
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D([—a,0]). First, we establish the Feller property for the Markov semigroup af- 
ter time t = a. Under the main assumption of a stable drift, we establish the 
tightness of the solution segments using semimartingale characteristics and apply 
the Krylov-Bogoliubov method to obtain an invariant measure on the Skorokhod 
space. Due to the absence of the strong Feller property Doob's method fails to 
prove uniqueness of the invariant measure. From an abstract point of view the 
loss of the strong Feller property is due to the degeneracy of the diffusion term 
when the equation is lifted to the segment space: the driving process is only one- 
dimensional, cf. Gatarek and Goldys [12] for the abstract non-degenerate case. 
The question of uniqueness of the stationary solution turns out to be subtle and 
the degeneracy of the noise process does not permit a straight-forward analytical 
treatment. While for certain cases uniqueness will be shown to hold, a counterex- 
ample lets us suspect that uniqueness fails in greater generality. Nevertheless, the 
correlation structure of the solution process, if it exists, is uniquely determined 
and analytically tractable. 

In the next section we briefly review some basic facts about stochastic delay dif- 
ferential equations. Section 3 is devoted to the variation of constants formula and 
properties of the Markov semigroup. The existence and uniqueness of stationary 
solutions are discussed in Sections 4 and 5, respectively. 

2 Preliminaries 

We follow standard notation, in particular we write C[a, b] for the space of real- 
valued continuous functions on [a, b]. The Skorokhod space of all real- valued 
functions on [a, b] that are right-continuous and have left limits at every point 
(cddldg for short) is denoted by D[a, b]. It is endowed with the Skorokhod metric 
d s given by 



where A[a, b] := {A : [a, b] — ► [a, b] : A is an increasing homeomorphism}. Note 
that ds(<p n ,<p) implies the convergences (p n {a) —* <p n (b) — ► p(b), but 

not the pointwise convergence in the interior (a, b). 

The space (D[a,b],ds) is a separable metric space. Moreover, there exists an 
equivalent metric d on D[a, b] such that (D[a, b], d) is a complete separable metric 
space, see for instance Jacod and Shiryaev [15]. We endow D[a, b] with the 
corresponding Borel cr-algebra B(D[a, b]). For ip € D[a,b] we denote by <p(t—) 
its left-hand limit at t and we define A</?(i) := (p(t) — tp(t—), t £ (a,b], and 
A<p(a) = 0. For a > and a function (p G D[— a, oo) we introduce the segment 
of ip at time t ^ as the function 



ip t : [—a, 0] — > R, tp t (u) := tp(t + u). 
Let us first turn our attention to the deterministic delay equation underlying the 
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stochastic equation (1.1): 



x(t) = <p(0) + f If x(s + u) n{du) ) ds fort^O, 

J \J[~a,0] J 



(2.1) 



x(u) = (p(u) for u E [— a, 0], 



where \i is a signed finite Borel measure and the initial function ip is in D[— a, 0]. 
Note that the inner integral exists because ip and a fortiori also x are measurable 
and locally bounded. 

As the fundamental system in linear ordinary differential equations and the Green 
function in partial differential equations, the so-called fundamental solution or 
resolvent plays a major role in the analysis of (2.1). It is the function r : R — > 1R 
which satisfies (2.1) with the initial condition r(0) = 1 and r(u) = for u E 
[— oo, 0). The solution x(-, ip) of (2.1) for an arbitrary initial segment p E D [— a, 0] 
exists, is unique, and can be represented as 

x(t, ip) = p(0)r(t) + f [ r(t + s-u)p(u)dufi(ds) for t ^ 0, (2.2) 
J[-a,0] Js 

cf. Chapter I in Diekmann et al. [11]. The fundamental solution converges for 
t — > oo to zero if and only if 

vo(jj) := sup J Re(A) : A E <D, A - / e Xs fi(ds) = I < 0, (2.3) 

{ Jl-afl] J 

where Re(z) denotes the real part of a complex number z. In this case the decay 
is exponentially fast and the zero solution of (2.1) is uniformly asymptotically 
stable. 

Let us fix a complete probability space (f2, P) with a filtration (Ftjt^-O satisfy- 
ing the usual conditions. We study the following stochastic differential equation 
with time delay: 



dX(t) = I f X(t + s) fi(ds) ) dt + F(X)(t-) dL(t) for t ^ 0, 

\J[-a,0] J 



(2.4) 



X(u) = $(u) for n E [—ot, 0], 



where ju is a signed finite Borel measure and the initial process ($(u) : n E 
[— a, 0]) is assumed to have trajectories in D[— a, 0] and to be jF -measurable. 
The driving process L = (L(t) : t ^ 0) is a Levy process. We denote its Levy- 
Khintchine characteristic by (b,a 2 ,v) with respect to the truncation function 

Turning to the specification of the nonlinear mapping F, we remark that results 
for the existence and uniqueness of strong or weak solutions of stochastic delay 
differential equations driven by Brownian motion appear in different generalities: 
Mohammed [21] provides a result under random functional Lipschitz conditions, 
Mao [18] discusses in addition the method of steps, which provides a unique 
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solution without a regular dependence of the coefficients on values in the past, 
Liptser and Shiryaev [17] give general results for weak solutions and ltd and Nisio 
[14] consider the existence of weak solutions for equations with finite and infinite 
delay. Since our equations are driven by Levy processes and the most general 
conditions are not our concern here, we follow Protter [24] and merely assume that 
the deterministic functional F : D[— a, oo) — > D[— a, oo) is functional Lipschitz 
and autonomous, i.e. it is continuous with respect to the Skorohod topology and 
it satisfies for all (pi G D[—a, oo), i = 1,2: 

(a) there exists a constant K > 0, independent of q>i and t, such that 

|F(¥>i)(t) - F(<f>2)(t)\ < K sup l^iCa) - p 2 (s)| for all f > 0; (2.5) 

(b) F(¥>i(a + ■))(*) = F{ip{)(t + s) for all t, s > 0. 

Equivalently, setting ^(^[-c^o]) := -^(vXO) the two conditions can be stated as 
F((p)(t) = F(ip t ) with a functional F which is Lipschitz continuous on D[— a, 0] 
equipped with the supremum norm. 

We can rewrite the differential equation (2.4) as the integral equation 

X(t) = &(0)+ [ G< s >(X)(s)ds+ [ H*(X)(s-)dL(s) for t ^ 0, (2.6) 
Jo Jo 

when introducing G^,, : D[0,oo) — ► D[0, oo) for s ^ and V G D[0, oo) by 
(abusing notation slightly) 

G*>(V0(«) = / M«) l[-a,- a )(«) + W«) l[- S) o](«)) /*(<*«)> (2-7) 

J[-a,0] 

i^(V)(*) = F {ip l[_ a)0 ) +V l[o,oo)) («). (2.8) 

For .^-measurable initial segments $ the mappings G<j> and H$ are functional 
Lipschitz in the definition of Protter [24] and we can invoke Theorem V.7 in 
Protter [24] which ensures a unique strong solution of (2.4). Recall that a strong 
solution of (2.4) is an adapted, stochastic process X with cadlag paths satisfying 
(2.6). The solution is called unique if all solutions are indistinguishable. We 
denote the solution by (X(i) : t ^ —a) or (X(t, <fr) : t ^ —a). 

Examples 2.1. 

(a) The no-delay case: if /j, = b8o, a point-mass at zero, and F(ip)(t) = f(f(i)), 
t ^ —a, then the equation reads 

dX(t) = bX(t) dt + f(X(t)) dL(t) for t > 0. 

If / is Lipschitz continuous, then F is easily seen to be functional Lipschitz 
and autonomous. 
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(b) The point-delay case: suppose fx = X)i=i an d F((p)(t) = f(<p(t — 
cti),...,ip(t - a n )), t ^ 0, and F(tp)(u) = F(ip)(0), u G [-a,0], with 
ai G [— a,0]. Then the equation reads 



n 



dX(t) = ^2biX(t - dt + f(X(t - ai), . . . , X(t - a n )) dL{t) for t ^ 
i=i 

and F is again autonomous and functional Lipschitz if / is Lipschitz in all 
its arguments. 

(c) The distributed-delay case: for p(ds) = b(s)ds and F((p)(t) = 

/(/[-a,0j¥>(* + a ) c ( a ) ds )> t > 0, and F(<p)(u) = F(<p)(0), u G [-a,0], we 
obtain for t ^ 

dX(t)= [ X(t + s)b(s)dsdt + f I [ X(t + s)c(s) ds) dL(t). 

J[-a,0] V^[" Q '°] / 

Again, we need / to be Lipschitz in order to have F functional Lipschitz 
and autonomous. 

(d) Further examples and counterexamples: other useful path-dependent map- 
pings like F(ip)(t) = sup Mg [ t „ Q)t ] f(u) and their combinations with Lips- 
chitz functions are functional Lipschitz and autonomous. Beware, however, 
that not all Lipschitz continuous functionals F on D[— a,0] give rise to a 
functional F : D[— a, oo) — > D[—a,oo), for instance the jump size func- 
tional F(tp)(t) = Aip(t) is not cadlag for cadlag functions ip with jumps. 
It is interesting to note that all admissible linear functionals are given by 
F(<p) = Jj_ a j <p(u) p(du) with p ranging through the space of finite Borel 
measures, which follows from the result by Pestman [23] when excluding 
the part based on jump sizes. 

3 Properties of the solution 

3.1 The variation of constants formula 

Many of our considerations will be based on a stochastic convolution equation, 
the variation of constants formula. This formula is easily derived if the driving 
process has bounded second moments, but no longer for processes where an ltd 
isometry or inequality fails. We provide a proof separately in Reifi et al. [26]. 

Theorem 3.1. Let F be functional Lipschitz. Then for a stochastic process 
X = (X(t) : t ^ —a) and initial condition <I> the following are equivalent: 

1) X is the unique solution of (2.4) with Xq = 

2) X obeys the variation of constants formula: 

x{t) = $ ) + I r (t ~ s ) F (X)(s-) dL(s), t > 0, (31) 

W), t€[-a,0], 
where r is the fundamental solution of equation (2.1). 
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3.2 Measurability of the segment process 



Our further work will be strongly based on considering the segment process (Xf : 
t ^ 0) in D[—a,0] instead of the real-valued process (X(t) : t ^ —a). This 
approach is natural because the segment process is Markovian and turns out to 
be eventually Feller. These properties will pave the way for our further analysis. 
In the case that L is a Brownian motion the segment process is immediately a 
Feller process on the path space C([— a, 0]), see Theorem III. 3.1 in Mohammed 
[21], which is not true in our setting because of the discontinuity of the shift 
semigroup on D[— a, 0]. 

The following two Lemmas establish certain measurability and continuity proper- 
ties of the segment process. For a continuous path space similar properties have 
been studied in Chapter 3.7 of Da Prato and Zabczyk [9] and Lemma II. 2.1 in 
Mohammed [21] in Chapter 3.7. 

Lemma 3.2. LetY = (Y(t) : t € [a, b]) be a progressively measurable, real-valued 
stochastic process with cadlag paths. Then Y is a D[a,b]-valued random variable. 

Proof. According to [5, Thm. 14.5] the Borel a-algebra of D[a,b] coincides with 
the cylindrical cr-algebra generated by all point evaluations ip *— > <p(c) for ip £ 
D[a, b] and arbitrary c 6 [a, b]. Hence Y is measurable as Y(c) is Borel measurable 
for every c £ [a, b]. □ 

Lemma 3.3. Let (Y(t) : t ^ 0) be a stochastically continuous process with 
cadlag paths. Then the segment process (Y t : t ^ a) in D[— a, 0] is stochastically 
continuous as well. Moreover, there exists a jointly measurable modification of 



Proof. For h > we define the homeomorphism : [—a, 0] — > [—a, 0] by A/ l (s) := 
s — h for s £ [—a + 2h, —h] and affine respectively on [—a, —a + 2h] such that 
Xh{— a) = —a and on [—h, 0] such that A^(0) = 0. Then \\Xh — Id ||oo ^ h and 



ds(Y t+h ,Y t ) < \\Y t+h o\ h - YfW^ + \\X h - Id ||oo -»• |Ay(i)| as h | 0. 



Hence, i i— > is right-continuous at to if Y is continuous at to- Similarly, one 
can establish lim/j|orfs(Yj_/ l , Yj) ^ \AY(t — a)\. We conclude that t ^ Yt is 
stochastically continuous at t if P(AY(t - a) ^ 0) = P(AY"(t ) / 0) = 0, 
which follows from the stochastic continuity of Y. 

Any stochastically continuous process with values in a Polish space has a jointly 
measurable modification, which is proved following [9, Prop. 3.2], but measuring 
the distance with the metric of this space. This gives the final assertion. □ 



(Yf.tZ a). 




s e [-a + 2h, -h], 

s E [-a, -a + 2h) U (-/t,0]. 



Therefore, we obtain 
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3.3 The Feller property 

Basic tools for deriving the existence and uniqueness of invariant measures are the 
Feller and strong Feller property of the Markov semigroup defined by the segment 
process. We establish here the Markov property, the Feller property after time 
a and give examples that the immediate Feller and the eventually strong Feller 
property fail in general. For our purposes the ordinary Markov property of the 
segment process will be sufficient, but the strong Markov property can also be 
derived following the lines of Chapter 9.2 in Da Prato and Zabczyk [9]. 

Proposition 3.4. Let X be the unique solution of (2.4). Then the segment 
process (X t : t ^ 0) is a Markov process on D[—a,0]: 

P(X t G B | T s ) = P(X t eB\X s ) P-a.s. 

for all t^s^O and Borel sets B £ B(D[-a,0]). 

Proof. We fix u ^ and consider for t ^ u the equation 

X u ^{t) = <p(0) + f I X u ^{s + v) n{dv) ds + f F(X u ^)(s-) dL(s), 

Ju J[—a,0] Ju 

X u ' lfi (m) = (p(m — u) for m E [u — a,u] and 92 € D[—a,0]. 

We denote the unique strong solution by (X u,ip (t) : t ^ u — a) and the segment 
process by (X^ ,ip : t ^ u). 

We define Q u := o~(L(s) — L(u) : s ^ u) which is independent of the cr-algebra 
T u from the given filtration. The solution X u -> Lp (t) is ^-measurable for every 
t ^ u and, as in the proof of Lemma 3.2, it follows that the segment X^ is 
^-measurable as well. The uniqueness of the solution implies X(s) = X u < Xu (s) 
for every s ^ u — a and thus Xt = X^' Xu for every t ^ u with probability one. 
By construction (cf. [24]) the solution process depends in a measurable way on 
the initial condition so that the function 

A : D[-a,0] x fi -> R, A((p,u) := l B (X?' v (w)) 

is measurable for every B € B(D[— a, 0]) and independent of T u for fixed (p. An 
application of the factorisation lemma [10, Prop. 1.12] yields P-almost surely 

P(x t eB\F u ) = ni B (x^ Xu ) \t u ] = nMXu,-) I T u \ = EL%v)]| v =x u , 

which ends the proof because the right-hand side is cr(X ?t )-measurable. □ 

Let Bb(D[— a, 0]) denote the space of all real-valued bounded Borel functions on 
D[— a, 0] endowed with the supremum norm IHI^ and Cb(D[—a,0]) its subspace 
of continuous functions. Due to Proposition 3.4 the operators 

P 8 ,t : B b (D[-a, 0]) - B b (D[-a, 0]), P s , t f( V ) := E[/(X t s '<0] 

have the property that P u ,sPs,t = Pu,t f° r ^ u ^ s ^ t. By homogeneity we 
have P Sj t = Po,t-s for ^ s ^ t, cf. Thm. V.32 in Protter [24], and the operators 
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P t : = Po,t) t ^ 0, form a Markovian semigroup. The Markovian semigroup will be 
called eventually Feller if there exists a to ^ such that for any / <G Cb(D[— a, 0]) 
the following two conditions are satisfied: 

P t f € C b (D[-a, 0]) for every t ^ t , (3.2) 
\imP s f{ip) = P t f(ip) for every <p € D[-a,0], t^t . (3.3) 

By definition the solution (X(t) : t ^ 0) has cadlag paths and it is easily ob- 
served from (2.6) that it has no fixed times of discontinuity: P(AX(t) ^ 0) ^ 
P(AL(t) / 0) = 0. Hence, the process (X(t) : t ^ 0) is stochastically continuous. 
By Lemma 3.3 so is the segment process (Xt : t ^ a) and thus condition (3.3) is 
fulfilled for to = ol. The semigroup is not stochastically continuous for to < a and 
condition (3.2) fails for to < ol due to the discontinuity of the shift semigroup, as 
the following example demonstrates. 

Consider < (3 < a and the initial functions ip n := l[- / 3(i- n - 1 ),o] which for 
n — ► oo converge in D[— a,0] to (p°° := Ij-^o]- The corresponding solution 
segments X™, for an arbitrary specification of F and L in the differential equa- 
tion, satisfy X™_J—a) = ip n (—f3) = 0, while X^_J— a) = 1 holds. Hence, 
ds(X™_p, X^_a) ^ 1, which implies that ip \— > Ptf(ip) is not continuous for 
/(^>) := a) | A 1 £ Cb(D[— a, 0]) and any time t € (0, a). Similarly, 
t i ^ Ptf(ip°°) is seen to be discontinuous at t = a — (3. 

We now establish condition (3.2) for to = ol by showing even more, namely 
that ip i — ^ (X(t, <^>) : t € [0, T]) is continuous from D[— a, 0] to the space of 
cadlag processes with the uniform convergence on [0, T] in probability, which is 
stronger than convergence in the Skorokhod topology in law. We start with a 
norm estimate in spirit of Emery's inequality before proving the main result. In 
accordance with Section V.2 in Protter [24] we employ the following norms for 
semimartingales (Z(t) : t ^ 0) and adapted cadlag processes (Y(t) : t ^ 0): 

\\Y\\ 2 s2[T] :=E[ sup Y(t) 2 ], 

\\Z\\ 2 H*[o,T] ■= inf {E[M, M] T + E[TV(,4)(T) 2 ]} , 

where the infimum is taken over all possible decompositions Z = M + A where M 
is a local martingale and A a bounded variation process with M(0) = A(0) = 0. 
The total variation of A on [0, T] is denoted by TV(A)(T). 

The quadratic variation process is defined by [Z,Z] := Z 2 — f Z(s-) dZ{s). 
Based on these norms the spaces H 2 [0,T] and S' 2 [0,T] are constructed canon- 
ically. Moreover, they are Banach spaces, and H 2 [0,T] is continuously embedded 
mS 2 [0,T}. 

Lemma 3.5. Suppose the Levy process L has a finite second moment and (H(t) : 
^ t ^ T) is an adapted cadlag process with j E[#(t-) 2 ] dt < oo. Then 

2 

H(s-)dL(s) 

H 2 [0,T] 



a 2 + Jx 2 is(dx) + (EI(l)) 2 rj| E[H(t-) 2 ]dt. 
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Proof. This follows from the decomposition L(t) = M{t) +t!EL(l) with M a 
square integrable martingale. □ 

The surprising result of the next proposition, which says that convergence of 
the initial conditions in Skorokhod metric implies uniform convergence of the 
solution processes, is essentially due to the fact that the driving Levy process is 
a semimartingale without fixed time of discontinuity. 

Proposition 3.6. Assume F : D[—a, oo) — > D[—a, oo) is continuous with respect 
to the Skorokhod metric. Let X n be the solution of equation (2.4) with determin- 
istic initial segment ip n , let (p n — > ip in D[— a,0] and let X be the solution with 
initial segment tp. Then (X n (t) : t ^ 0) converges to (X(t) : t ^ 0) uniformly on 
compact sets in probability. 

Proof. We consider first a stopping time R such that the process L R ~ is a- 

sliceable for some suitably small a > in the sense of [24]. 

In analogy to [24, Thm. V.10] we use the representation (2.6) and put 

Y n {t) := j\G v {X)-G ip n{X)){s)ds 
Jo 

+ f (H v (X)-H^(X))(s-)dL R -(s), 
Jo 

g n {U){t) := G^(X){t) - G^(X - U)(t), 
H n {U){t) := H v n(X)(t) - H v n(X - U)(t), 

t ^ 0, to obtain for U n := X — X n the equation 
U n (t) = 99(0) - v? n (0) +Y n (t) + f g n (U n ){s)ds+ f H n {U n ){s-)dL R -{s). 



By [24, Lemma V.3.2], extended to two driving semimartingales, the solution 
U n of this equation satisfies ||t/ n || 52 [ t] ^ ^ llvW ~~ ^"(0) + ^ n |ls 2 [oT] f° r an y 
T > with a constant C > depending on the process L R ~ and a uniform 
bound for the Lipschitz constants of G^n and H^n . The Skorokhod metric ensures 
V? n (0) -» (p(0), so that ||^ n || 5 2 [0jT] -> follows if Y n tends to zero in S 2 [0,T]. 
The latter is fulfilled if 

e[ / (G v (X)(t)-G^(X)(t)) 2 + (H ip (X)(t)-H ip n(X)(t)) 2 dt\ (3.4) 
Jo 

tends to as n — > 00, due to the continuous embedding H 2 S* 2 
from [24, Thm.V.2] and Lemma 3.5 with the additional observation that 
\\J Q J(s-)dL R -(s)\\ H 2 [0tT] s$ \\J Q J(s-)dL(s) ||tf2[o,T] for any process J e5 2 [0,T], 
Let w be fixed for the moment. The functions H v n(X(uj)) converge in the Sko- 
rokhod topology to H v (X(uj)), which implies convergence in L 2 [0,T]. Concerning 
the sequence G Vn we have 

[ T (G Vn (X(u,))(t)-G ip (X(u,))(t)) 2 dt= [ [ f T (<p n (t + u) 

JO J[—a,0) J[—a,0) Jo 

<f(t + u))(<p n (t + v ) - ip(t + v)) l[_ a - t )(u) l[_ a _ t )(u) dtfi(du) fi(dv). 
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This expression converges to zero as n — > oo, since the Skorohod convergence 
of ip n to p implies (p n — > p Lebesgue a.e. and the sequence (<^ n )n is uniformly 
bounded. Again by [24, Lemma V.3.2] the solution process A is an element of 
S 2 [0,T], whence by the uniform linear growth of (G^™) and (H^n) the argument 
inside the expectation in (3.4) is dominated by a P-integrable function. The 
Dominated Convergence Theorem thus gives the convergence in (3.4) such that 
\\X - X n \\ s2[0yT] = ||£/ n || 52[0iT] -► for any T > 0. 

Next, let L be an arbitrary Levy process. According to [24, Theorem V.5, p. 192] 

there exist stopping times = To ^ T\ ^ T2 ^ • • • such that sup^ Ti = 00 a.s. and 
L Te- 

is a-sliceable for each I. Consider equation (2.4) with L replaced by L Te 
and let X n ' e denote the solution with initial segment p n and let A°°' £ denote 
the solution with intial segment ip, for n, £ G N. We have shown above that 
X n ^ — > X 00 ' 1 uniformly on compact sets in probability for every I. Further, it is 
clear from the equation that X n ' e = (X n ) Te ~ . Let now t > 0, r > 0, and e > 
be arbitrary. Choose an I such that P(Te < t) < e/2. Then 

P( sup |A n (s) - A(s)| > r) 

< P( sup \X n {s) - X(s)\ ^ r and T £ ^ t) + P(T e < t) 

s$ P( sup \X n ' e {s) - X°°' e (s)\ >r)+e/2<e 

for n large. Hence X n — > A uniformly on compact sets in probability. □ 

Let us finally show that in general we cannot expect that the solution is eventually 
strongly Feller, which is characterised by the existence of a to > such that for 
all f e B b (D[-a,0}) 

P t f G C b (D[-a, 0]) for every t ^ t . 

Using indicator functions for /, this implies 

p 1 ^ P(X t (tp) G P) G C b (D[-a, 0]) for every B G B(D[-a, 0]), t ^ t . 

Suppose the functional F in the equation (2.4) is of the form F(ip)(t) = f(ip(t—a)) 
for t ^ and F(tp)(t) = for t < 0, with a Lipschitz-continuous homeomorphism 
/ : R — > (a, 6), 6 > a > 0, and consider the case that L is standard Brownian 
motion. Then the quadratic variation (X t ) of the solution segment X t , t ^ a, 
satisfies 

/•t+u 

(X t ) u = / f 2 (X (s - a)) ds P-a.s. for u G [-a, 0]. 
Jo 

Since both sides of the equation are continuous in u for continuous A, there is 
one P-null exception set for all u G [—a, 0]. Consider the map 

^)H : =/ _1 ((^r) 1/2 )' «e[-a,0], 
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defined on the functions tp with finite quadratic variation such that ^ G (a 2 , b 2 ) 
for Lebesgue-almost every u G [— a, 0]. We have P(ypQ)(u) = X(t + u — a), u € 
[— a, 0]) = 1 for all t ^ a. Iterating this map, we can recover with probability 
one the initial segment Xq from observing X ma since V m {X ma ) = Xq for every 
integer m. This identifiability property shows that the laws of the segments 
X ma ((fi) and X ma (if2) for different initial segments ipi and <£>2 must be singular. 
Hence, there is a contradiction to the strong Feller property at to = which 
asserts the continuous dependence of the laws on the initial condition. 
In fact, this example even shows that the Markov semigroup is not eventually 
regular in the sense of Da Prato and Zabczyk [10]. We shall see in Section 5.2 
that this counterexample is due to the delay in the diffusion coefficient. 

4 Existence of a stationary solution 
4.1 Tightness 

We establish the tightness of the laws {C(X t )}t^o i n D[0, a] by considering the 
semimartingale characteristics. Recall that (6, a 2 , v) denotes the Levy-Khintchine 
characteristic of the Levy process L. 

Assumption 4.1. 

(a) The delay measure (i in the drift satisfies vo(fi) < with vq from equation 
(2.3). 

(b) The jump measure v satisfies JLi^ log \x\ v(x) < oo. 

(c) The coefficient F in equation (2.4) is functional Lipschitz, uniformly 
bounded and autonomous. 

Condition (a) yields the exponential decay of the fundamental solution, while 
condition (b) ensures that J * f(s) dL(s), for exponentially decaying functions 
/ of locally bounded variation, converges in law and is already for constant F 
necessary for the existence of a stationary solution, as was shown by Gushchin and 
Kiichler [13], cf. also Thm. 4.3.17 in Applebaum [2]. In condition (c) restrictions 
on F are imposed such that the differential equation is autonomous, has a unique 
solution and the impact of the driving process cannot become too large. For the 
latter the imposed boundedness of F can certainly be relaxed considerably, but 
will then depend on the large jumps of L, that is, on fine properties of v. 

Proposition 4.2. Grant Assumption 4-1- Then the solution process (X(t) : t ^ 
—a) of (2.4) with initial condition Xo = has one- dimensional marginal laws 
{£(X(t))}t^o that are tight. 

Proof. Let us split the Levy-process L into two parts, one of them consisting of 
jumps of size larger than one: 

L{t) = N(t) + R{t) with N(t) = AL(a) 1{|al( s )|>i} • 
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Then the variation of constants formula (3.1) yields X = Y + Z with 

Y(t):= [ r(t- s)F(X)(s-)dN(s), t ^ 0, 
Jo 

and Z(t):= [ r(t - s)F(X)(s-) dR(s), i > 0. 
Jo 



Tightness of (X(t) : t ^ 0) will follow from tightness of Y and Z. 
The fundamental solution r decays exponentially with \r(t)\ ^ ce _/3i for some 
constants c, (3 > due to Assumption 4.1(a). Considering Y first, we obtain for 
any K > with m := sup^, \F(tp)(0)\ by time reversal for the compound Poisson 
process N the estimate 



P(\Y(t)\ > K) < P |r(t - s)F(X)(s-)AiV(s)| > A 
< P \J2 ce ~ m ~ s)m \ AN ( s )\ > K 



s<t 



K 



= P Ee-*|AJV(,)|> — 

The tightness for Y follows from the tightness of X^<t e_/3s |AA(s)| which has 

been established in [13, Lemma 4.3] under Assumption 4.1(b). 

Since R is a Levy process with bounded jumps its canonical decomposition is, by 

means of [2, p. 103], of the simple form R(t) = R (t) +FER(1) where (Ro(t) : 

t ^ 0) is a square-integrable martingale. We split Z into the sum Z = Zq + Z\ 

with 

Z {t):= f r{t-s)F{X)(s-)dR {s), t> 0, 

Jo 

and Zi(t) := E[i?(l)] / r(t - s)F(X)(s-) ds, t ^ 0. 

For Zi we easily obtain P(|Zi(t)| > K) ^ P{cmf3- 1 | E[P(1)] | > if) = for A" 
sufficiently large, implying tightness of Zi. As in Lemma 3.5 we obtain 

E[Z (t)] 2 < ( a 2 + f x 2 u(dx))m 2 f r 2 (t-s)ds, t ^ 0. 

\ V|x|<l / JO 

Hence by the exponential decay of r, the sequence (Zo(i))t^o is bounded in Lp(f2) 
and thus tight. □ 

Proposition 4.3. in £/ie setting of Proposition we have that the laws 
{C(X(t + s) - X(t), s € [0,a])} m are tight in D[0,a\. 
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Proof. We are led to consider for t ^ and s £ [0,a] 



Y t (s) :=X(t + s)-X(t) 

rt+s 



^ du + J t + F(X)(u-)dL(u). 



X(u + v) n(dv) 

-a,0] 

Let us introduce for t ^ the semimartingale (It(s) : s € [0, a]) by letting 

rt+s 

h(s):=J F(X)(u-)dL(u) for a €[0, a]. 

Now, either by following the lines in [15, III.2.c] and using the Levy-Ito de- 
composition or by applying [25, Prop. 7.6] the semimartingale characteristic 
(Bj t ,Ci t ,ui t ) of (It(s) : s € [0, a]) is found to be 

rt+s 

+ J xF(X)(«-)(l ( _ lil) (xF(X)(«-))-l ( _ 1)1) (x))i/(da;))d« ) 

/•t+s 

C It (s)=a 2 J F 2 (X)(u-)du, 
vj t (ds, dx) = ds x Kj t (X, t + s, dx), 
where for y G £>[— a, oo), u ^ t and a Borel set A € £>(IR) 

K It (y,u,A) := J l A \{o} (F(y)(u-)x) v(dx). 

Hence, the semimartingale (Y t (s) : s € [0, a]) has the characteristic (By t , Cy t , vy t ) 
with Cy t = Ci t , vy t = vj t and 



B Yt (s) = B It (s) + [ t+S If X(u + v) n{d- 

Jt \J[-a,0] 



v) du. 



We prove the tightness of (Y t ) t ^ by means of [15, Thm. VI.4.18] and [15, VI.4.20]. 
For that we have to verify that 

a Yt (s) :=TV(B Yt ){s) + CY t (s)+ f (\x\ 2 A 1) vy t (du, dx), s€[0,a], 

forms a tight sequence (ay t )t^o of processes and all limit points of the sequence 
{£(ay t )}i^o as t — > oo are laws of continuous processes. According to [15, 
Prop. VI. 3. 33 and VI. 3. 35] this will follow if there exist some increasing pro- 
cesses (Ay t )t^o satisfying these conditions and in addition Ay t — ay t defines for 
every t ^ an increasing process, since ay t (s) ^ a.s. for all s G [0, a], t ^ 0. To 
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obtain such processes (Ay t (s) : s € [0, a]), we estimate 



TV(B Yt )(s) 



t+s 



■a,0] 



X(u + v) n{dv) + bF(X)(u-) 



+ / XF{X){U-) (!(_!,!) K^) 



-a,0] 



+ \b\ m + c du, 



a 



where m := sirpy, \F(ip)(0)\ and the finite constant c is defined by 
c := / m |x| i/(dar) + z/(R\(-l, 1)). 

Vi<b|<l 

Therefore, the process ay t is majorized by 

A Y t {s):= ! [ X(u + v)n(dv) 
Jt J[-a,0] 

+ s(\b\m + c + (T 2 m 2 + J ((m 2 x 2 ) A 1) u(dx)^j for s £ [0, 

and ^4y t — ay ( is increasing. Since -Ay t is continuous and only depends on t in the 
first term, it suffices to prove tightness in C[0, a] of the first term: 

rt+s r 

Jy t (s) := / / X(u + v) /j,(dv) du for s E [0, a]. 

J[-a,o] 

Recalling r(u) = for u < 0, we obtain by the variation of constants formula 

J(u) := / X(u + u) /i(du) 
■/[-a,o] 

= / (I r(u + v-s)F(X)(s-)L(ds)] fi(dv) 

J[-a,0] \Jo J 

f'U 

= / r(u- s)F(X)(s-)dL(s). 
Jo 



To prove tightness of the absolutely continuous processes (Jy t )t^o> it suffices to 
show that the process (I (t+s) : s G [0, a])t^o is bounded in probability in C[0, a], 
that is, 



lim supP ( sup \I(u)\ > K ) =0. 

K^oo t j> \t<^u^t+a 



(4.1) 



Note that we have the exponential decay estimate \r(t)\ ^ c'e _/3t . Decomposing 
L into its drift, diffusion, and large and small jump parts, it is clear that only 
integration with respect to the large jump part N may pose problems. As in the 
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proof of Proposition 4.2, however, the restriction on the large jumps in L and the 
finite intensity of N yield in a similar manner 



sup 



/ r(u- s)F(X)(s-)dN(s) < ^ cW - ^* - *) \AN(s 



and the tightness of the right-hand side by [13, Lemma 4.3]. Thus we infer 
(4.1). □ 

Theorem 4.4. Grant Assumption 4-1- Then for the solution process (X(t) : t ^ 
—a) of (2.4) with initial condition Xq = the laws of the segments {C(X t )}t^ a 
are tight in D[— a,0]. 

Proof. If we let Z t (s) := X(t — a) for s £ [— a, 0], then the processes (Z t )t^o of 
constant functions are tight in C[— a, 0] by Proposition 4.2. On the other hand, 
{C(Xt — Zt)}t^ a are tight in D[— a, 0] by Proposition 4.3 applying the time shift 
t i — ► t — a. Therefore the sum (X t — Z t ) + Z t is tight in D[— a, 0] using the result 
in [15, VI.3.33(a)]. □ 

4.2 From tight solutions to stationary solutions 

We use the construction due to Krylov and Bogoliubov, see for example Da Prato 
and Zabczyk [10]. For the reader's convenience we include a complete proof, 
which is tailored for our purposes. Consider equation (2.4) and its Markovian 
semigroup (Pt)t^o as defined below Proposition 3.4. Denote by V = V(D[—a,0]) 
the set of Borel probability measures on D[—a,0], endowed with the topology 
of weak convergence of measures. Let (•, •) denote the duality pairing of V and 
B b := B b (D[-a,0}) given by ((J) = J fd(, (eV,feB b . Define for t ^ and 
C € V the functional Pj*( by 

(Pt*0f:=(C,Ptf), f€B b . 

If ( is the distribution of an initial segment <£, then P t *£ is the distribution of 
X t ($), since 

(P t *C,/> = J E[/(X t (^))]C(^) =E[E[/(X t ($))|^ ]] =E[/(X t (d>))], 

for / € B b . A measure ( £ V is called an invariant measure or stationary 
distribution of (2.4) if P t *( = ( for all t ^ 0, that is, ((,P t f) = ((J) for all 
f e B b and all t ^ 0. 

It follows from Lemma 3.3 that t >— ► is a continuous map from [a, oo) to V 
and moreover P* +t ( = PgPtC f° r s,t ^ 0. Further, P t maps C b := C b (D[-a,0}) 
into C b for all t ^ a, by Proposition 3.6. 

Because of Theorem 4.4, the next theorem follows from Theorem 4.6 below. 

Theorem 4.5. Grant Assumption 4-1- Then there exists a stationary distribu- 
tion for (2.4). 
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Theorem 4.6. If for some £ G V the set {-P t *C ■ t ^ a} is tight, then there exists 
an rj G V such that P^n = n for all t ^ 0. Moreover, n is an element of the closed 
convex hull of {PfC ■ t ^ a} inV . 

Proof. Denote for convenience T t := Pt+ a and C{t) := T^n, t ^ 0. First we show 
that for each t ^ there exists a unique i5 f 6? such that 

<<?*,/> = 7 f\as),f)d8 for all fe C b . 
T Jo 

It is routine to show the uniqueness. In order to show existence, define 

V>(f) ■= T /'(COO, /> ds for every / G C b . 
1 Jo 

Then ip G , where C£ denotes the dual Banach space of C b . Let e > and take 
# C D[-a,0] compact such that ((s)(H) ^ 1 - e for all s € [0,t]. For / G C fe 
with H/lloo ^ 1 and /(x) = for all x G i7 we then have 



KCW,/>I 



<ll/lloo C(s)P[-a,0]\ff)<e 



for all s G [0, t], so | <£>(/) | ^ e. Now by the Riesz-Bourbaki representation theorem 
(see for instance [7, Prop. 5.2.5] and [7, Prop. 5.6.12]), there exists a tight finite 
positive Borel measure i? on D[— a,0] such that <p(f) = for every / G C b . 

Notice that (■&, l D [- afi ]} = 1, so i? S P. 

Next we show that i? is an element of the closure of the convex hull of {C(s) : ^ 
s ^ t} in "P. Let M denote the weak* closure of the convex hull of {C(s) : ^ 
s ^ i] in C b . Then M is a weak* closed convex set in C b . The Hahn-Banach 
Theorem implies that for any tp G C£ \ M there exist f £ Cb and /? G R such 
that < 13 and (r?, /) > /? for all 77 G M. Then (0, /) = \ jJ<C(s), f) ds > 0. 

Thus ■& G M and therefore iJeMnP. 

Since {C(s) : s ^ 0} is tight, its convex hull is tight and hence relatively compact 
in V by Prohorov's Theorem. Thus the set {dt '■ t ^ 0} is contained in a compact 
set and therefore there exist a sequence t n f 00 and a measure n £ V such that 

#t„ -»• »?■ 

Finally, for £ ^ a and / G C& we have 

1 /"in 1 rtn 

lim - / {C(t + s),f)ds= lim - / (T:,P t f)ds=( V ,P t f) 
and on the other hand 

i r n ( C (t+ S ),/)d S = l [ tn ((( S )j)ds-±- f\as)j)d8 

In JO r n JO r " JO 

1 rtn+t 

+ -/ (COO, /><**, 



which converges to (77, /) as n — ► 00. Hence, P t *r/ = 77 and it follows that P^rj = 
Pt(P*rj) = Pt +a n = rj, for every t ^ 0. □ 
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We remark that the proof given above remains true in a more general setting. 
Indeed, we need only to replace D[— a, 0] by an arbitrary separable metric space 
E and assume that (Tf)t^o is a family of bounded linear operators on Cb(E) 
such that T s +t+a = T s Tt for all s,t ^ 0, that for some ( £ one has that 

T^Q € V(E) for all t ^ 0, that the map t \— > (T t *£, /) is measurable from [0, oo) 
to R for all / G C 6 (£), and that the set {TfC : t ^ 0} is tight. 

5 Uniqueness of the stationary solution 

As we have seen, the Markovian semigroup is in general not eventually strongly 
Feller so that a main tool to establish uniqueness of the invariant measure is 
not available. Moreover, when considered as a stochastic evolution equation, the 
generator of the deterministic equation (2.1) is only eventually compact (see [11]) 
and the Markov semigroup is only weakly continuous with a generator which is 
analytically not easily tractable (see [21]). Hence, typical analytical methods to 
prove uniqueness (see [19] for a survey) cannot be easily applied either. 
We therefore consider several specific cases where uniqueness can be proved nev- 
ertheless: for small Lipschitz constants by a contraction argument, in the Wiener 
case for non-delayed diffusion coefficients by establishing the strong Feller prop- 
erty via Girsanov's theorem and for compound Poisson driving processes and 
non-delayed drift terms by studying the deterministic behaviour between the 
jumps. After that, we relax the requirements and show that in full generality 
second-order uniqueness holds up to a constant factor. We conclude by an exam- 
ple where the invariant measures are not unique. 

5.1 Small Lipschitz constants 

If the function F is not too far from being constant, as measured by the Lipschitz 
constant, then uniqueness holds. The upper bound for the Lipschitz constant 
below can be reconstructed by our proof, but it is certainly not the best possible. 

Theorem 5.1. Grant Assumption 4-1 and suppose that the Levy process has 
finite second moments. If the Lipschitz- constant K of F in (2.5) is sufficiently 
small then the laws of all stationary solutions X of (2.4) coincide. 

Proof. Let X and Y be two stationary solutions with corresponding initial con- 
ditions Xo and Yq. As mentioned above Proposition 5.7 the moments IE ||-Xo|Ito 
and E || lo|| are finite. 

As vo([i) < the fundamental solution r decays exponentially with |r(£)| ^ ce~^*, 
Jo°° r2 ( s ) exp(2/3s) ds < oo, and J °° r 2 (s) exp(2/3s) ds < oo for some constants 
c, (3 > 0. We choose an arbitrary constant 7 < (5 and we let Z{v) := F(X)(u) — 
F(Y)(u) for convenience. By use of the decomposition L(t) = M(t) + E[L(1)] t 
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with a martingale M, the variation of constants formula implies for t ^ a: 



E 



sup \e^ s (X(s) -Y(s))\' 

jt— a^.s^.t 



< 3E 



sup |e 7S x(s,Xo-y )r 

.t— a^s^t 



sup 



+ 3E 

t— a^J 

+ 3(E[L(1)]) 2 E 



e 7S r(s - u)Z(u-) dM{u) 

! e^ s r{s-u)Z{u-) 
Jo 



(5.1) 



sup 

t—aSjs^t 



du 



An application of representation (2.2) yields 



E 



sup \e^ s x{s,X -Y )\' 

jt— a^.s^.t 



^d®\\X -Y \\' c 



(5.2) 



for a finite constant d depending only on the measure fi. Let r\ be the function 
defined by r\(s) := r(s) exp(7s). Then we obtain for the second term in (5.1) 



E 



sup 



f e^ s r{s-u)Z{u-)dM{i 
Jo 









E 


sup 






t— a^s^t 





< 2E 



sup 

t— a<s<t 



n (0) + / h (m) dm ) e< u Z(u-) dM{u) 
\ Jo 



e~< u Z{u-) dM{u) 



+ 2E 



sup 

t—a<s^t 



'0 V0 

The first term in (5.3) can be estimated by 



e lu Z(u-)dM{u) him) dm 



(5.3) 



E 



sup 



e lu Z{u-)dM(u) 
a 2 + J x 2 is{dx)^ j\^ u -E\Z(u-)\ 2 du. 



(5.4) 



Note that f\ has essentially the same asymptotic as s h exp(~fs)r(s). Hence 
if we choose a constant 6 > such that 7 + <5 ^ we obtain d\ : = 
J °° exp(25m) |ri(m)| 2 dm < 00. Applying Holder's inequality to the second term 



19 



in (5.3) results in 



E 



sup 

t— a<s<t 



ai; 



< E 



sup f 

t-a^s^t JO 



e< u Z(u-) dM{u) him) dm 



|fi(m)| dm 



/ / e^ u Z(u-) dM{u) 
Jo Jo 



-25m 



dm 



^ d ie - 26 ^ E 



sup / / 

-a^s^tJO JO 



e^ u Z(u-) dM(u) 



e 2Sm dm 



< d ie 25a (a 2 + j x 2 v(dx)^j e- 25m drn^j J* e 2 ^ E \Z(u-)\ 2 du. (5.5) 

The last term in (5.1) can be estimated similarly by Holder's inequality 



E 



sup 

t—a<s<t 



I e ys r(s — u)Z(u—) du 
Jo 

^°°e 27U |r(n)| 2 duj jf* e 2 ' yu E \Z(u-)\ 2 du. (5.6) 



By collecting the inequalities (5.2) to (5.6), using the Lipschitz condition (2.5) 

1 1 2 

and applying Gronwall's Lemma we conclude E \\X S — Y s \\ — > for s — > oo if 
the Lipschitz constant K is sufficiently small. Consequently, the laws of Xo and 
Yq coincide. □ 

5.2 Non-delayed diffusion coefficient 

We have seen in Section 3.3 that the Markov semigroup [Pt)t^o of the solution 
segments is in general not eventually strong Feller. This is only an effect due 
to the delay in the diffusion term and cannot be caused by a delayed drift for 
the Wiener-driven case, as we shall see now. Let us consider as special case of 
equation (2.4) 

dX(t)=\[ X(t + s)u(ds)) dt + f{X{t))dW{t) fort^O, (5.7) 

\J[-a,0] J 

with initial segment $ as in (2.4), a Wiener process W and a Lipschitz function 
/ : R — > R. By a simple argument based on Girsanov's theorem we obtain the 
following result. 

Proposition 5.2. If f satisfies the ellipticity condition mi x ^ f(x) > 0, then 
the solution segments (X t : t ^ 0) of (5.7) generate a Markov semigroup on 
C([— o,0]) that is strongly Feller after time a. 

Proof. First note that the continuous functions form a closed subspace of the Sko- 
rokhod space D[— a, 0] such that the formerly obtained results are in the Wiener- 
driven case also valid on C[— a, 0]. Referring to Theorem 7.19 for diffusion-type 
processes in [17], we infer from the Lipschitz continuity of the coefficients and 
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from the ellipticity of / that the laws Q\ and Q2 of the solution processes of 
(5.7) on C[0, T], T > arbitrary, are equivalent for different delay measures /ii 
and fj.2- The corresponding Radon-Nikodym derivative is given by 

|g(X) = ex P (^ T (| a0] X(t + s) f(X(t))- 2 dX(t) 

o] ^(* + a)(Mi-M2)(<fe)j /(X(t))- 2 dt). 

As in [20, Thm. 2.1] one can show that the validity of the strong Feller property 
at each time is invariant under the change of measure. According to that result 
we need to check that the semigroup is Feller and that 



km E 



dQ \x\^))- d S±{x\.^)) 



dQi dQi 



= 



for initial segments <p n — > <p in C[—a, 0] and for the corresponding solution process 
X 1 with the choice //i. The Feller property has been established in Section 3.3. 
By Scheffe's Lemma it suffices for the second condition to prove convergence in 
probability. This is accomplished by the continuity of the map tp 1— ► X l {-, ip) from 
C[-a,0] to L 2 ([0,T] x 0) for any T, which follows from [21, Thm. 3.1]. 
We have thus reduced the problem to proving the strong Feller property of the 
Markov semigroup generated by the solution segments (X t )t^o of 

dX(t) = f(X(t)) dW{t) for t ^ 0, (5.8) 

as special case of (5.7) with jjl = 0. It is well known that this diffusion equation 
generates a strongly Feller semigroup on R under our assumptions on /, see e.g. 
[10, Thm. 7.1.1]. We claim that this property is inherited by the segment process. 
For this consider a bounded measurable functional VP on C[—a, 0] and remark that 
X(-;<p) = X(-;<p{0)) only depends on the initial value, not the whole segment. 
By the scalar Markov and weak uniqueness property we obtain for t ^ a and any 
initial segment <p with obvious notation 

E[tt(X t (y>))] =nn^(Xt^))\J r t-a}} = ^[^[^{X a {X{t - a-^u),J))]). 

Setting H{£) := E[\l/(X a (£))], £ € R, the scalar strong Feller property implies 
the continuity of 

r) 1 ^ B[H(X(t - a; V ))} = E w [E w ,[*(X Q (X(t - a; r?,u;),u/))]] 

for 77 G R. Since (p n — ► ip in C[— a,0] yields ip n (0) — > y?(0), we thus infer the 
continuity of 

¥> -> E w [E w /[*(X Q (X(t - a;ip,uj)),uj')]} = E[*(X t {<p))] 
on C[— a, 0], which is the asserted strong Feller property at t ^ a. □ 
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Corollary 5.3. The Markov semigroup (Pt)t^o ^ regular after time 2a. Thus, 
any stationary solution of (5.7) is unique and strongly mixing. 

Proof. Recall that we have regularity at to if all transition probabilities 
P(X to (ip) G •) are equivalent for ip G C[— a, 0]. By Doob's Theorem [10, Thm. 
4.2.1] this property yields the uniqueness and strong mixing result. 
The regularity property at to > 2a is implied by the strong Feller property 
at a together with the irreducibility at to — a [10, Prop. 4.1.1], which means 
that all transition probabilities at time to — a have support in the entire space. 
To prove the latter, we may again restrict to the case ji = and consider X 
as in (5.8) due to the equivalence of the laws. As in [27, Cor. VIII. 2. 3] it 
follows from Girsanov's theorem that the support of the (regular) conditional 
law £(X to _ a \X(t - 2a) = x) is given by S x := {/ G C[-a,0] : f{-a) = x}. 
Since the law of X (to — 2a; <p) has for the same reasons the full support R for any 
initial segment (p, we conclude by composition that £(X to - a (<p)) has full support 
C[— a, 0] independent of tp, which yields the required irreducibility. □ 

5.3 Uniqueness in the compound Poisson case 

Let us consider here the case of a Levy triplet (b,a 2 ,u) with a = 0, b = 0, and 
the total variation A := H^H^y finite, that is L is a compound Poisson process. 
If there is no delay in the drift, then we can reduce the question of uniqueness 
of the invariant law on the Skorohod space D[—a,0] to a property of the one- 
dimensional invariant law. 

Proposition 5.4. Suppose L is a compound Poisson process and consider for 
a > a differential equation of the form 

dX{t) = -aX(t)dt + F{X)(t-)dL{t) fort^O, (5.9) 

admitting a strong solution for any initial segment. If an invariant solution mea- 
sure on D[— a,0] exists and the one- dimensional marginal distributions of any 
two invariant measures are non-singular, then the invariant measure is unique. 

Proof. Let pi and p2 be two invariant measures. By coupling methods we can 
construct a filtered probability space (O,, P, (Pt),P) carrying the process L, and 
the .Fo-nieasurable random variables Y,Ze D [— a, 0] with P Y = pi, P z = p2 and 
P(Y(0) = Z(0)) > 0. Denote by X 1 and X 2 the corresponding strong solution 
processes with initial conditions Y and Z, respectively. 

Since with probability e~ Aa > the process L does not jump on the interval 
[0, a], we have 

p(xi = xl y(o) = z(o)) ^ p(£ ! A ^)I = o,y(o) = z(o)) 

= P(£|AL(t)|=0)P(y(0) = Z(0))>0. 

Hence, introducing the set 

S := {cp\3u; G n : X^(u) = X 2 a {u) = <p} C D[-a,0], 



22 



we find for any Borel set B in D[— a, 0] 

P(X l a £BnS)> P({u | X 2 a (co) = Xfa), Xl(u) G B}) 
= P({u>\Xl(u)=X 1 a (u;),Xl(u)eB}) 
=■ r s (B) 

and equivalents P(X% e B S) ^ t s (B). By invar iance, we conclude 

mm{p 1 (Br\S),p 2 (BnS)} ^ t s {B) for all B G B(D[-a,0]) 

with a non-negative measure T5 satisfying ts(S) > 0. Hence pi and p2 are non- 
singular; for if pi(A) = and p 2 (A ) = for some Borel set A, then 

rs(S) = t s (S nA) + t s (S n A c ) <^ Pl (snA) + P2 (s n A c ) = o. 

As extremal points of the set of invariant measures are singular (see [10, Prop. 
3.2.7]), uniqueness follows. □ 

Theorem 5.5. Grant Assumption 4-1 ■ Suppose L is a compound Poisson process 
and consider equation (5.9) with a > and F ((f)(0) > for all f G D[—a,oo). 
Then there exists a unique invariant measure for (5.9). 

Proof. If the jump measure v is zero, then the compound Poisson process vanishes 
and the only invariant measure is clearly the point measure in zero. Let us now 
first consider the case of possible positive jumps: ^((0, 00)) > 0. By Proposition 
5.4 it suffices to show that any two invariant one-dimensional distributions are 
non-singular. We first show that they are absolutely continuous with respect to 
Lebesgue measure. 

Let B CM, denote any Borel set. For the solution process X of (5.9) we find 
P(X(t) G e- at B) > P(J2 |AL(a)| = 0, X(0)e~ at G e~ at B} 

= e~ xt P(X(0) G B). (5.10) 



Now assuming that X is stationary with one-dimensional marginal law po, we 
obtain by Fubini's Theorem for any Lebesgue null set B and T > 

/ po(e~ at B)dt= I I t e -at B (x) p (dx) dt 
Jo Jo Jr 

= / y l B (t)dtp o (dx) + Po ({0})l B (0) 

JR\{0} Jx 1 

= p o ({0})t B (0). 

By estimate (5.10), however, the left-hand side is bounded from below by 
po(B)^^ — . Hence, we infer po(B) = for all Lebesgue null sets B with ^ B. 
Since F is positive and v 7^ 0, we can exclude a point mass in zero because the 
state {0} will be eventually left by the process P-a.s. and the probability to jump 
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back exactly to this state is zero. We conclude that po is absolutely continuous 
with respect to the Lebesgue measure. 

Let S denote the support of po. Since F is positive and bounded away from zero 
and L has positive jumps, there will occur with positive probability sufficiently 
many positive jumps of L in short time that the trajectory X will take arbitrarily 
high values. This means for the support S of the marginal invariant measure po 
that sup S = +00. 
For a Borel set B C (0, 00) we have 



If Po (B) = 0, then (5.10) with B replaced by e at B yields that p (e at B) ^ 
e xt po(B) = for all t > 0, and we obtain that 



Since sup S = +00, we infer that the Lebesgue measure of B equals 0. Thus the 
Lebesgue measure of (0, 00) is absolutely continuous with respect to p$. 
If v((— oo,0)) is also positive, then the symmetric argument yields that the 
Lebesgue measure on R is equivalent with po. In any case, we know that two 
invariant measures are both equivalent to the appropriate Lebesgue measure and 
hence with each other. An application of Proposition 5.4 completes the proof. □ 

Remark 5.6. We have derived the regularity property that, unless the jump mea- 
sure is zero, the one- dimensional marginals are absolutely continuous with respect 
to the Lebesgue measure. 

In some cases one can easily derive the density of the invariant measure. For 
example, if we assume L to have only positive jumps of size at least J > and 
F ((f)(0) £ [<7o,0"i] for all (f G D[— a, 00) and some uq, o\ > then the density of 
the marginal of the invariant measure of (5.9) near zero is given by 



with a suitable constant C . 

5.4 Second-order uniqueness 

A real-valued stochastic process (X(t) : t ^ —a) will be called second-order 
stationary, if < E[A(i)] 2 < 00, the values E[X(i)] are constant for all t ^ 
—a, and the function (s,t) 1— > E[X(s)X(i)] depends only on the difference s — 
t. Obviously, any stationary solution of (2.4) with finite second moments is 
second-order stationary. If the Levy process is a square-integrable martingale, we 
establish second-order uniqueness for equation (2.4) up to a constant factor, more 





Six) = C-x 



(\—a)/a 



X e [0, Jcr ), 



a 
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precisely the expectation and the correlation function are uniquely determined 
and can be calculated analytically. 

Note that the invariant measure exhibited in Section 4 will have finite second 
moments for its one-dimensional marginal whenever the Levy process has finite 
second moments. This follows from the fact that the constructed tight sequence of 
segments (X t ) will be uniformly bounded in L 2 P {Q) by the variation of constants 
formula (3.1) and Lemma 3.5. 

Proposition 5.7. Grant Assumption 4-1- Suppose the Levy process is a square- 
integrable martingale with characteristics (b,a 2 ,v). Then any stationary solution 
(X(t) : t ^ — a) of (2.4) with finite second moments is a centered random process 
with auto-covariance function 

c(h) := lE[X(0)X(h)} = Var p (0)] r(s)r{s + h) ds, h ^ 0. 

H r lli2(R + ) J0 

The spectral density is given by 

e » E[A(0) 2 ] (||r|| i2(R+) \xM)\y 2 > ^R. 

where Xn( z ) := z ~ I[- a o] e zu ^{du) is the characteristic function of the determin- 
istic equation (2.1). 

Proof. By the variation of constants formula (3.1) and the martingale property 
of L we have for t ^ 

EI(t) = Ej(t,I o ) = r(i)E[I(0)]+ / f r(t + s - u)E,[X(u)]du n(ds). 

J[-a,0] Js 

Due to limj^oo r(t) = and stationarity we conclude that X is centered. Again 
using the variation of constants formula, we find for h, t ^ 

ft 



7E{X(t)X{t + h)]=m[x(t + h,X ) f r{t-u)F{X){u-)dL{u)} 

Jo 

rt+h 

+ E,[x(t,X ) / r(t + h-u)F(X)(u-)dL(u)] +E[x(t,X ) x(t + h,X 
Jo 

rt rt+h 

+ E[/ r(t-u)F(X)(u-)dL(u) / r(t + h - u)F(X)(u-) dL(u)\. 
Jo Jo 



As in Lemma 3.5 we obtain 

r t rt+h 



rt rt+n 

/ r(t-u)F(X)(u-)dL(u) / r(t + h - u)F(X)(u-) dL{u) 
Jo Jo 

a 2 + J x 2 v(dx) S j r(t-u)r(t + h-u)^[F{X)(u-)} 2 du. 



E 



The variance is estimated as the expectation before: 
Var[x(i,X )] 



2 ( r(t) 2 Var[A(0)] + ( / f \r(t + s - u)\ E \X(u)\ du \n\(ds) 

\ \J[-a,0] J-s , 



25 



which converges to as t — > oo. Applying the Cauchy-Schwarz inequality to the 
first three terms in the equation above results in 



Cov(X(0),X(h)) = lim TE[X(t)X(t + h)] 



t— »oo 



E[F(X)(0)] 2 (a 2 + j x 2 u(dx)^j jf r(u)r(u + h) du. 



This yields the expression for the covariance function. The formula for the 
spectral density follows from the fact that r is the inverse Fourier transform 
of Xn(—i-)~ l , as obtained for affine stochastic delay differential equations driven 
by a Wiener process in [16]. □ 

Remark 5.8. It is seen from the proof that 

VarLY(O)] = E[F(X)(0)] 2 [a 2 + | x 2 v{dx)\ \\r\\ 2 L2{R+) , 

which gives some information about the size of the variance depending on bounds 
for the functional F. We shall see in the counterexample of Section 5.5 that this 
variance term need not be uniquely determined, at least for measurable functionals 
F. 



5.5 Non-uniqueness 

In the Wiener-driven case we construct an elliptic diffusion functional F which 
remains constant in time for certain initial segments, but with different values for 
different initial segments. By doing so, we can recover, for instance, Ornstein- 
Uhlenbeck processes with different diffusion coefficients as solutions. Suppose F 
is of the form 

F{<p){t) := ^max(l, min(|(^)*:^ /2 , 2)) t u+ (t) for t > -a, 

where {ip) b a denotes the quadratic variation of tp € D[— a, oo) on the interval [a, b] 
which might be infinite. Then F is bounded away from zero and infinity and 
is measurable (as a limit of measurable functionals), but obviously not continu- 
ous. Leaving slightly our framework, let us consider for a Wiener process W the 
equation 

dX{t) = -X{t)dt + F{X)(t-)dW{t) fort^O, 
X{u) = $(u) for u G [-a,0]. ^ ' ' 

Due to the positive minimal delay a/2 there exists a strong unique solution 
to this equation for any .Fo-measurable initial segment by the method of steps, 
cf. Mao [18]. On the other hand, there exists for every a £ [1,2] a stationary 
Ornstein-Uhlenbeck process X a which solves the equation (we suppose that W 
is a two-sided Wiener process) 

dX a (t) = -X a {t)dt + adW{t) for i€lR. 
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Then choosing $> a = Xq , we obtain that each X a is also a stationary solution 
of (5.11). This is due to the fact that (X CT }*I° /2 = § a 2 and thus F(X a )(t) = a 
hold for all t ^ and a £ [1, 2]. 

This example shows that some kind of regularity of F has to be imposed to 
guarantee uniqueness, but we do not know whether already for functionals F with 
large, but finite Lipschitz constants uniqueness breaks down. It is interesting to 
note that a similar dichotomy has been described by Mohammed and Scheutzow 
[22] for the long time behaviour in dependence of the diffusion functional. 
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